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ABSTRACT 

We have examined the radial distribution of white dwarfs over a single 
HST/ACS field in the nearby globular cluster NGC 6397. In relaxed popula- 
tions, such as in a globular cluster, stellar velocity dispersion, and hence radial 
distribution, is directly dependent on stellar masses. The progenitors of very 
young cluster white dwarfs had a mass of ~ 0.8 M©, while the white dwarfs 
themselves have a mass of ~ 0.5 M . We thus expect young white dwarfs to 
have a concentrated radial distribution (like that of their progenitors) that be- 
comes more extended over several relaxation times to mimic that of ~ 0.5 M & 
main-sequence stars. However, we observe young white dwarfs to have a signif- 
icantly extended radial distribution compared to both the most massive main 
sequence stars in the cluster and also to old white dwarfs. 
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Introduction 



The large proper motions of some pulsars indicate that they have space velocities of 
hundreds of kilometers per second, which were presumably imparted to them when they 
became neutron stars (for recent reviews see iPodsiadlowski et al.ll2005l ; lRomanill2005l ). The 
question naturally arises, could something of this sort have happened to white dwarfs? They 
are, of course, not observed to have comparably high velocities, and their births do not 
involve anything as energetic as a supernova, but there is already a hi nt in the literature 
that white dw a rfs ma y have higher velocities than their progenitors did. IWeidemannl (119771 ) 
and IWilliamsl (120021 ) have observed that the numbers of white dwarfs observed in open 
clusters are lower than expected — a natural explanation being that a fracti on of the white 



dwarf s have escaped from the cluster. To investigate this question further, iFellhauer et al. 



( 120031 ) performed N-body simulations and found that a kick of magnitude approximately 
twice that of the cluster velocity dispersion (several kms" 1 ) would deplete a loosely bound 
open cluster of almost all its white dwarfs. 

The hypothesis that white dwarfs begin their lives with a small velocity excess leads 
to scenarios that are particularly conducive, for dynamical reasons, to testing in a globular 
cluster. The progenitors of the newest white dwarfs in a typical globular cluster had a mass 
of ~ 0.8 M Q , but they have become white dwarfs with a mass of ~ 0.5 M . Assuming 
a quiescent birth, they begin their lives with the spatial distribution that is appropriate 
for their progenitors, but over the course of a relaxation time — of the order of a few times 
10 8 years — they will acquire the more extended spatial distribution that goes with their 
new, lower mass. However, if white dwarfs are given natal kicks, they then start their lives 
with a velocity dispersion that is too large for their original spatial distribution, and in the 
course of a crossing time — only a million years or so — they acquire a more extended spatial 
distribution, perhaps even more extended than they will have after relaxation has made 
their spatial distribution correspond to their new masses. In a typical globular cluster, stars 
well down the white dwarf cooling sequence have been white dwarfs for dozens of relaxation 
times, while the youngest white dwarfs have entered that stage very recently. Moreover, 
these clusters have hundreds of white dwarfs, so that the test can be made in a statistically 
significant way. Finally, because globular clusters are old stellar populations, the mass 
difference between the progenitors of the youngest and oldest white dwarfs in our sample is 
< 0.09 Mm- Though the white dwarf initial-final mass relation is uncertain (IFerrario et al. 



20051 ; iKalirai et al.l 120051 ) . the mass difference between the white dwarfs must certainly be 
much less than the mass difference between their progenitors. This allows us to treat all 
white dwarfs as essentially of equal mass. 



The question then is, do young white dwarfs in a globular cluster have a spatial dis- 



- 3- 



tribution that is more extended than that of their progenitors? Until recently, the problem 
with addressing this question has been that white dwarfs in a globular cluster, particularly 
the old, fainter ones, are at a magnitude where they cannot be easily observed. Today, 
however, deep imaging with HST is able to re ach faint enough to image the entire white 
dwarf sequence in very nearby globular clusters ([Richer et al.ll2006l ; lHansen et al.ll2007l ). We 
will show here that in NGC 6397 the radial distribution of its young white dwa rfs is signifi- 



cantly more extended than that of their progenitors and also old white dwarfs. iDavis et al. 



(2006]) presented preliminary results on the white dwarf distribution in this cluster and (jHeyl 



2007a| jbl) then used Monte Carlo simulations to generate radial distributions of stars with 



kicks at birth and found that a kick of the order of 1.8 tim es the veloci t y disp ersion of the 
cluster giants was required to account for the effect seen by IDavis et al.l (120061 ) . 



Observations 



In NGC 6397 we used HST's ACS Wide Field Channel to image a single field centere d 
5' SE of the cluster center for 126 orbits (GO-10424, for details see IRicher et al.l 120061 ). 
The exposure time was divided between the F814W and F606W filters. The field ranges 
from 179" to 391" from the cluster center, corresponding to 1.3-2.8 half-mass radii ((Harris 
19961 ). Proper motions, for the elimination of field stars, were measured on archival WFPC2 
images. We were able to extend our magnitude limit considerably by blind measurement 
of the WFPC2 images at every point where the ACS images told us that there was a star, 
but the shorter exposure time of the WFPC2 images was still what set the limit of our 
completeness. Even so, artificial-star tests showed that at the magnitude that corresponds 
to a white dwarf age of 5 Gyr our completeness was still 92%. Table [TJ summarizes the major 
properties of NGC 6397. 

The essence of this Letter is to compare the radial distributions of white dwarfs of various 
ages with those of main-sequence stars of various masses. Neither mass nor age is directly 
observable; in what follows we will need models of both main-sequence stars and cooling 
w hite dwarfs. To assi gn a mass to each main-sequence star, we use main-sequence models 
of iBaraffe et al.l (J1997J), with the metallicity appropriate for NGC 6397. To assign an age to 
each white d warf, we use a set of c ooling models that wer e developed using interiors from 
Wood! (119951 ) and atmospheres from iBergeron et al.l (119951 ) . Our cooling models are for 0.5 
M & stars with pure carbon cores and thick hydrogen-rich atmospheres. Shown in Figure [TJ 
is the colour-magnitude diagram (CMD) of NGC 6397 with the main-sequence masses and 
white dwarf cooling ages superimposed. 
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Radial distributions 



To test whether two distributions differ significantly, we make use of two statistical tests: 
the familiar Kolmogorov-Smirnov (KS) test, which uses the maximum difference between the 
normalized cumulative distributions, and the Wilcoxon rank-sum (RS) test (equivalent to 
the Mann- Whitney U-test), which compares rankings within the merged, sorted distribution. 
Each of the tests gives the probability that two random samplings from the same parent 
population will show differences as large as the observed ones, though the RS test (because 
it can be applied as a one-tailed test) is generally more sensitive. 



Mass segregation is known to be present in NGC 6397 (IKing et al.lll995t lAndreuzzi et al. 



20041 ) . We test the sensitivity of our data set to mass segregation by comparing the radial dis- 
tributions of 0.15-0.2 M & main-sequence stars with 0.5-0.6 M main-sequence stars (i.e. the 
presumed true mass of the white dwarfs). Figured left section, shows the cumulative radial 
distributions of these two groups, and the power with which the KS and RS tests can dis- 
tinguish between them. Our data are clearly able to detect mass segregation, even over the 
limited radial extent of our field. 

The white dwarfs that have had time to undergo dynamical relaxation should have a 
distribution similar to that of the ~ 0.5 M & main-sequence stars, while the unrelaxed ones 
should have a more centrally concentrated distribution, like their more massive progenitors. 
However, NGC 6397 exists in the Galactic potential, and periodic changes in the lumpy 
potential could affect the re laxation process. In fact, the orbit of NGC 6397 has a disk- 
crossing time of ~ 100 Myr (IMilone et al.ll2006l ; iKalirai et al.l 120071 ). and one might expect 
disk shocking to play a role in the internal cluster dynamics. However, in relaxed clusters, 
such as NGC 6397, we expect these effec ts to be most impo rtant close to the tidal radius, 
and insignificant at the half-mass radius (IGnedin et al.lll999l ). 



To study the distributions of the white dwarfs, we need a "young" and an "old" sample. 
(When we refer to age here, we mean the time since the star became a white dwarf.) We 
want all the stars of the young sample to be young enough not to have experienced appre- 
ciable relaxation since they acquired the lower mass that goes with being a white dwarf. 
Therefore, we need to look at the relaxation time in our region — bearing in mind, however, 
that re l axatio n time is only an approximate time scale rather than an exact quantity. The 
Harris! (119961 ) catalog gives a relaxation time of 0.29 Gyr at a half-mass radius of 2'.3 for 
NGC 6397. Given that even the inner edge of our field is at a radius where the star density 
is lower than it is at 2f3, our stars should remain unrelaxed for ~ 1 Gyr. 



For choice of the old-sample boundaries there are two considerations. The lower bound- 
ary in age should be high enough so that the stars have had enough time to have relaxed to 
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the distribution that is appropriate for their white dwarf mass. The upper age limit must be 
low enough to ensure that the completeness at the magnitude of the faintest white dwarf is 
not a function of crowding, and therefore radial position. We take the old group to be from 
1.4-3.5 Gyr (62 stars). 

For the choice of the young-sample boundaries there is a different set of considerations. 
The lower boundary should not be less than a crossing time. However, this timescale is of 
the order of Myrs — shorter than the youngest inferred white dwarf age in our sample — so 
this boundary can be set to zero Gyr. In principle it would be advantageous to set the upper 
boundary to select only the very youngest (i.e. age much less than a relaxation time) white 
dwarfs for this sample, thereby ensuring that the young white dwarfs would have had very 
little time to undergo interactions with other stars. In practice, it is necessary to set the 
upper boundary to an age comparable to a relaxation time in order to increase the sample 
size, while maintaining a sample largely free from the effects of relaxation. We choose a 
young sample that goes from 0.0-0.8 Gyr (22 stars). We note that we are in no danger 
of inducing a false signal by choosing an upper boundary that includes some relaxed white 
dwarfs, but we will rather dilute any signal that is truly there. 

While our definitions of young and old were guided by physical principles and observa- 
tional constraints, they remained somewhat arbitrary. In fact, there is some sensitivity to the 
upper boundary of the young group. If that boundary is placed anywhere from 0.4 to 0.8 Gyr 
the result remains strongly significant, increasing the boundary beyond this decreases the 
significance. Figure [21 right section, compares the radial distributions of the young sample 
with the old one in NGC 6397. 

Contrary to what one would expect supposing a quiescent birth scenario, the younger 
white dwarfs tend to be farther from the cluster center than old white dwarfs. This suggests 
that they have a higher velocity dispersion than old white dwarfs. However, this is not the 
key point. The most direct test would be to compare the radial distribution of young white 
dwarfs with that of their progenitors. Obviously, this test is impossible, as main-sequence 
stars of the mass of the progenitors no longer exist in the cluster. Instead we will test the 
white dwarf distributions against a series of main-sequence mass samples. Assuming that 
all main-sequence stars are relaxed, and that the trend is monotonic with mass, if we can 
show that the radial distribution of the young white dwarfs is inconsistent with the radial 
distribution of the highest mass main-sequence stars that we observe, we have shown that 
young white dwarfs are not drawn from the same population as their progenitors. 

We use the RS test to calculate the probability that the radial distributions of the young 
and old white dwarfs could come from samplings of the same parent distribution as main- 
sequence stars of various masses; the results are shown in Figure [3j Setting the acceptance 
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level at 0.05 (2a in a one-tailed normal distribution), we can reject the hypothesis that the 
young white dwarfs have a radial distribution consistent with that of main- sequence stars 
more massive than ~ 0.5 M & in NGC 6397. Furthermore, we note that we cannot reject 
the hypotheses that these same stars have radial distributions consistent with that of lower 
mass main-sequence stars. The rejection of the former hypothesis leads us to conclude that 
the young white dwarfs in NGC 6397 have a radial distribution more extended than that 
of their progenitors. Lending confidence to our method, Figure [3] also shows that the radial 
distribution of the old white dwarfs is not inconsistent with that of main-sequence stars of 
their supposed mass, but is inconsistent with that of low mass (M < 0.25 M & ) main-sequence 
stars. 



4. Are white dwarfs born with a kick? 

We have shown that the radial distribution of young white dwarfs in NGC 6397 is 
more extended than would be expected given a quiescent birth. Other than a possible 
kick at birth, is there an alternate explanation for this observation? One might imagine 
a mechanism related to the white dwarfs being in binary systems, either through binary 
disruption or the inhibition of white dwarf formation through c ommon envelope evo lution 



or mergers. However, the binary fraction in NGC 6397 is low flCool fc Boftonll2002l ). and 
any explanation invoking binaries is implausible. Alternatively, the newly born white dwarfs 
could undergo very close encounters during their first passage through the inner regions of 
the cluster with their reduced masses, and have their velocities increased to greater-than- 
equilibrium values through stellar interactions alone. We have produced N-body simulations 
that accurately model both binary evolution and stellar encounters, but omit natal kicks; 
we do not observe the young white dwarfs to have extended radial distributions in thes e 



simulations. A full description of these N-body models can be found in lHurley et al.l (120071 ). 
NGC 6397 has a collapsed core, and one may consider if interactions with the central density 
cusp could somehow create this effect. However, the white dwarf progenitors experienced 
the same potential as white dwarfs, so it is difficult to see how this could have an effect. 

Assuming the existence of a natal kick, we can estimate its size in the following way. At 
the radial distances that we are c onsidering, the measur ed velocity dispersion of ~ 0.8 M 



giants in NGC 6397 is 3.3 kms 1 ( iPryor fc Meylanlll993l ). In accordance with equipartition 



of energy in a relaxed system, we expect the velocity dispersion of stars of mass M to be given 
by 

f"vel — OViant 

^/M~jM. Thus, cr ve i for 0.2 M Q stars should be about twice that of the 
giants. Because velocity dispersions add in quadrature, the magnitudes of impulsive kicks 
that would transform the radial distribution of a 0.8 M Q population to that of a 0.5 M & or 0.2 
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M & population are 3 kms 1 and 6 kms 1 , respectively. These values are comparable to the 
velocity dispersion of a typical globular cluster, but still well below their escape velocities. 

The fact that we are examining a relatively small range of radii, and that we have a 
small number of stars in the white dwarf sample, prevents us from putting tight constraints 
on the velocity dispersion of the young white dwarfs. However, we note that the young white 
dwarfs have an extended radial distribution, implying an extra source of velocity dispersion. 
The most natural explanation for this observation is the existence of a natal kick. We leave 
it to the theorists to search for a detailed mechanism, but note that asymmetric mass loss 
while the st ars are on the asympto tic giant branch stages could generate a kick of the inferred 



magnitude (IFellhauer et al.ll2003l ) 
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Table 1. NGC 6397: Major Properties 



Property 


Value 


Distance Modulus (true) 


12.03 ±0.06 1 


Af814W 


0.33 1 ' 2 


[Fe/H] 


-2.02 ±0.07 3 


R-half mass 


140" 4 


thalf mass 


0.29 Gyr 4 


Observed field radius 


179-391" 



References. 



(2) ISirianni et al.1 
fl2003h : f4) lHarris 



1) Hansen et all (120071): 
fl2005h : f3 N ) lKraft fc Ivans 
Jl996h 



- 11 - 



14 


i i 


i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | 


1 1 1 1 


16 




/ 0.8 M© 


_ 








18 




\. . 0.6 M 




20 




3^ 0.4 M© 








'it 0.2 M© 




22 








24 


/ 


0.1 M© 

• 0.5 Gyr " r —-,.. . 






\ 


.1.0 Gyr 




26 


- \ 


• 1 2.5 Gyr 

: . 6.0 Gyr 

9 ** . 




28 




' ' 1 1 1 1 


■ III 



1 2 3 

F606W-F814W 



Fig. 1. — The CMD of NGC 6397. The white dwarf ages, according to the Wood/Bergeron 
white dwarf cooling models, and the main- sequence masses, as determined from models of 
Baraffe, are indicated. The magnitude range of the old (lower) and young (upper) white 
dwarfs are shown with arrows. 
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Fig. 2. — Left section: A comparison of the radial distributions of 0.5-0.6 M and 0.15-0.20 
M Q main-sequence stars for NGC 6397. The upper panel shows the RS test and the lower 
panel the KS test. The vertical lines in the upper panel indicate the radial position of each 
star (0.5-0.6 M Q above and 0.15-0.20 M below). Right section: A comparison of the radial 
distributions of the young and old white dwarfs. The panels are the same as in the left figure 
with the vertical bars denoting the locations of the old (upper) and young (lower) white 
dwarfs. Surprisingly, the young white dwarf population is less centrally concentrated than 
the old population. 
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Fig. 3. — The RS probability that the radial distribution of a given white dwarf population 
could be drawn from the same parent distribution as the radial distribution of a given mass 
of main-sequence star for NGC 6397. The dashed lines indicate the equivalent of one-tailed 
1<7, 2<t, and 3<r results. Note: the mass value of the points are the mean values of the 
mass bins. Adjacent points making up the lines are not independent, and are composed of 
approximately 90% the same stars. 



